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Abstract 



The suppression of hadron production originated from the induced jet energy loss is one of 
the most accepted and well understood phenomena in heavy ion collisions, which indicates the 
formation of color deconfined matter consists of quarks, antiquarks and gluons. This phenomena 
has been seen at RHIC energies and now the first LIfC results display a very similar effect. In 
fact, the suppression is so close to each other at 200 AGeV and 2.76 ATeV, that it is interesting 
to investigate if such a suppression pattern can exist at all. We use the Gyulassy-Levai-Vitev 
description of induced jet energy loss combined with different nuclear shadowing functions and 
describe the experimental data. We claim that a consistent picture can be obtained for the 
produced hot matter with a weak nuclear shadowing. The interplay between nuclear shadowing 
and jet energy loss playes a crucial role in the understanding of the experimental data. 
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1. Introduction 

Energy loss of high energy quark and gluon jets penetrating dense deconfined matter 
produced in ultrarelativistic heavy ion collisions leads to jet quenching and thus probes 
the quark-giuon plasma formed in those reactions [1-8] . The non-abehan radiative energy 
loss suppresses the hadron production yield in the momentum range 2 — 3 GeV/c < px < 
15 — 20 GeV/c at RHIC energies [9-12]. Recent preliminary experimental data on charged 
hadron production collected in PbPb collisions at LHC energies, especially at 2.76 ATeV, 
display the same phenomena [13]. Comparing the nuclear modification factor, Raa, at 
RHIC and LHC energies, one can see that the measured values are very close to each other 
in the momentum range 5 GeV/c < pt < 15 GeV/c. This finding is very much surprising, 
since the measured charged hadron multiplicities display a factor of 2 enhancement at 
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LHC energy with respect to the values at RHIC energies [14]. Thus we expect a two 
times higher entropy production, which can be connected to a two times higher density 
of colored parton density formed in the central heavy ion collisions. 

In order to investigate the influence of radiative energy loss on hadron production, we 
apply a perturbative QCD (pQCD) based description of heavy ion collisions, including 
energy loss prior to hadronization. First, we check that the applied pQCD description 
reproduces data on charged production in p + p collision at LHC energies. Our results 
are based on a leading order (LO) pQCD analysis. Detailed discussion of the formalism 
is published in Refs. [15-17]. We performed next to leading order calculations [18], the 
results will be indicated in pp collisions, however, for simplicity we would like to discuss 
our results at LO and compare the experimental data. 

2. Parton model for pp collisions 

Our pQCD calculations follows the general description of charged hadron produc- 
tion [15-17]: 

ahcd 

Here we use LO parton distribution functions (PDF) from the MRST98 parameter- 
ization [19] and a LO set of fragmentation functions (FF) [20]. The applied scales are 
Q — kpt/ and Q ~ npx- the NLO calculations are performed on the basis of Refs. [18]. 
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Fig. 1. Charged hadron yield in pp-coUisions at ^/s = 7 TeV measured by the CMS Collaboration [21]. 
The calculation is performed in LO (full line) and NLO (dashed line), see details in the text. 
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Utilizing available high transverse- momentum (5 < pr < 100 GeV/c) p + p data on 
charged hadron production at 7 TeV [21] we can determine the best fitting description. 
Fig. 1 shows a general agreement between the data and our calculations in LO and in 
NLO. The proper scale choice is k = 1 in LO and k = 2 in NLO calculations. Interesting to 
note that we do not need to introduce intrinsic-fcT^ at LHC energies, which ingredient was 
an essential part of proper description of hadron production at RHIC energies [17,18,22]. 
In fact we keep the earlier extracted intrinsic-fcy value at RHIC energies, but neglect it 
at LHC energies. 



3. Nuclear effects in heavy ion collisions 



Now we turn to the calculation for AA collision at RHIC and LHC energies. We include 

the isospin asymmetry and the nuclear modification (shadowing) into the nuclear PDF. 
We consider the average nuclear dependence of the PDF, and apply a scale independent 
parameterization with the shadowing function Sa/A{x)'- 



fa/A{x) = Sa/A{x) 



^ fa/p{x) + (l - f ) fa/n{x) 



(2) 



There are different shadowing functions to be applied, such as the relative weak old 
HIJING-shadowing function [23] and its newer version with and without 6-dependence [24]. 
This new version contains a very large suppression at low-x, which can be checked by 
the new LHC results. Another version, namely the EKS-shadowing function has been 
developed with a Q^-dependence but with a weak suppression at low-x, especially in the 
high-Q^ domain. In this calculation we will use the EKS99 [25], which is the basis of the 
later version and it is a LO-set. 

At RHIC energy we follow our earlier calculations and include intrinsic-fcr and nuclear 
multiscattcring [17,22,27], however, at LHC energies we neglect these contributions. 

In our calculation jet quenching reduces the energy of the jet before fragmentation. 
We concentrate on j/cm = 0, where the jet transverse momentum before fragmentation is 
shifted by the energy loss, p*{L/\) = pc — AE{E, L). This shifts the Zc parameter in the 
integrand to z* = zd (1 — AE/pc) [26]. The applied scale in the FF is similarly modified, 
Q = kpt/z*, while for the elementary hard reaction the scale remains Q = npc- 

With these approximations the invariant cross section of hadron production in central 
A + A collision is given by [18] 

Eh^^ = J d% (fr tA{r) tA{\h - r|)x 

fdXidX2fa/A{xi,Q'^)fb/A{X2,Q'^) ^'^ ^h/c{Zc,Q ) 



'ibcd 



dt Zc TTZc 



Here tA{b) = J dz pA{b, z) denotes the nuclear thickness function, and it is normalized 
as usual: / d^btA{b) = A. In case of heavy nuclei the Wood-Saxon formula is applied for 
the nuclear density distribution, pA(b- z). The integral in b indicates the nuclear overlap 
in central collisions and the consideration of the Glauber geometry. 
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4. Averaged jet energy loss 



Let us summarize our basic knowledge about non-abelian energy loss in hot dense 
matter. First estimates [1,2] suggested a linear dependence on the plasma thickness, L, 
namely /S.E w 1 — 2 GcV(L/fm), as in abclian electrodynamics. In BDMS [3,4], how- 
ever, non-abelian (radiated gluon final state interaction) effects were shown to lead to 
a quadratic dependence on L with a larger magnitude of AS. Similar results have been 
obtained from light-cone path integral formalism [5,6]. 

In the GLV formalism, applying opacity series (see Refs. [7,8]), finite kinematic con- 
straints were found to reduce greatly the energy loss at moderate jet energies. On the 
other hand, the quadratic dependence on L has been recovered in wide energy region. 

Non-abelian energy loss in pQCD has been calculated analytically in two limits. In 
the "thick plasma" limit, the mean number of jet scatterings, n = L/X, is assumed to 
be much greater than one. For asymptotic jet energies the cikonal approximation applies 
and the resummed energy loss (ignoring kinematic constraints) reduces to the following 
simple form [3-6]: 

ATP - 93^ r, cx^ 
^^BDMS = — 7 — V , (o) 

where Cr is the color Casimir of the jet (= Nc for ghions), and f-i^/Xg oc is a transport 
coefficient of the medium proportional to the parton density, p. The factor, t; ~ 1 — 3 
depends logarithmically on L and the color Debye screening scale, /x. It is the radiated 
gluon mean free path, A^, that enters above. 

In the "thin plasma" approximation [7,8], the opacity expansion was applied and in 
the first order the following expression was derived for the energy loss: 



/ 





dV^e// 2k^-q^(k-q)iL2 



7r(qi-F/z2)2 16x2 i;2 + (k-q)4i2 



Here the opacity factor L/Xg is the average number of final state interactions that the 
radiated gluons suffer in the plasma. The upper transverse kinematic limit is 

k^,, = min [4E^x\4E^x{l - x)] , (5) 

and the upper kinematic bound on the momentum transfer is q^^^^ = s/4 ~ SEii. 

Furthermore, ii^jj/fi^ = 1 + M^/^max- -^t RHIC energies, these finite limits cannot be 
ignored [7,8]. The integral averages over a screened Yukawa interaction with scale /i. The 
integrand is for an exponential density profile, p oc exp{—2z/L), with the same mean 
thickness, L/2, as a uniform slab of plasma of width L. 

It was shown in Refs. [7,8] that in the asymptotic E oo limit, the first-order expres- 
sion (4) reduces to the BDMS result [3,4] up to a logarithmic factor \og{E/p). Numerical 
solutions revealed that second and third order in opacity corrections to eq. (4) remain 
small (< 20 %) in the kinematic range of interest. 
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Fig. 2. Calculation of nuclear suppression for charged hadrons at RHIC energy of y/200 AGeV, including 
different shadowing functions. Wc choose opacity L/X = 5 (left panel) and L/X = 6 (right panel). Exper- 
imental data are from PHENIX Collaboration for charged particle (square) [10] and tt^ (triangle) [11]. 



5. Numerical results 



We have calculated jet energy loss at RHIC and LHC energies and compared our 
simplified model description to existed data. We applied the jet quenching effects for a 
generic plasma with an average screening scale fi — 0.5 GeV, as = 0.3, and an average 
gluon mean free path Ag = 1 fm. 

In Fig. 2 we display our results for charged hadrons including three types of shadowing 
(HIJING, EKS99, new-HIJING and b-dependent new-HIJING) and jet energy loss at 
RHIC energy y/s = 200 AGeV. The opacity which reproduces the experimental data is 
between 5 and 6, we may choose L/X — 5.5. This result is consistent with the gluon 
rapidity densities of dNg/dy = 1000 — 1400 [28,29] and the charged hadron rapidity 
densities 550 — 650) [30,31]. One can see the weak dependence on nuclear shadowing, 
applying all three versions. 

Fig. 3 displays the results in PbPb collisions at ^/s = 2.76 ATeV, applying opacities 
L/X — 5 and 6. We can claim that using the new-HIJING shadowing [24] the requested 
opacity is also between L/X = 5 and 6, so L/X = 5.5 could give also a proper description. 

However, this result indicates that both at RHIC and LHC energies we have the same 
densities of the colored deconfined matter to be produced. This result contradicts the 
total charged multiplicity values measured at RHIC and LHC energies [14]. 

Thus we increased the opacity and assumed a larger parton densities. Fig. 4 displays 
results, where the ALICE data are reproduced at L/X — 7, but with the old HIJING 
and the EKS99 shadowing. This means a larger parton density because of the quadratic 
dependence on opacity. 

In fact the experimental data on Raa is so close together at RHIC and LHC, that 
there is no room for a strong shadowing, jet energy loss explains the suppression pattern. 
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Fig. 3. Calculation of nuclear suppression for charged hadrons at LHC energy of = 2.76 ATeV, 
including different shadowing functions. Wc choose opacity L/X = 5 (left panel) and L/A = 6 (right 
panel). Experimental data are from ALICE Collaboration for charged particle (dots) [13]. 
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Fig. 4. The nuclear modification factors at opacity value L/A = 7 at different shadowing functions, 
compared to ALICE data at LHC energy [13]. The good agreement with the old HIJING [23] and 
EKS99 shadowing [25] is emphasized. 
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Finally, we would like to emphasize the nice agreement between ALICE data [13] and 
the jet quenching scenario at high-p^, where the energy loss becomes weaker and suppres- 
sion factor approaching unity. The flat nuclear modification factor seen at RHIC energies 
disappear. That flatness was the results of a nice interplay between nuclear shadowing 
and multiscattering with the jet quenching. As the flat nuclear modification factor dis- 
appear at LHC energy, then the simple geometrical explanation of nuclear modification 
factor also becomes obsolete. This statement is one of the strongest consequences of the 
ALICE data. The other consequence is the negligible role of shadowing at LHC energies. 

Our analysis focused on induced jet energy loss and nuclear shadowing. The introduc- 
tion of elastic energy loss beyond the radiative energy loss may change the results, but 
the basic conchision on the presence of a weak nuclear shadowing will not be modified. 

Furthermore, these calculations clearly indicate the perturbative QCD region, where 
the parton model including jet energy loss is applicable to explain the measured hadron 
suppression pattern in heavy ion collisions. This region is > 6 — 7 GeV/c. Under this 
momentum value the behaviour of the nuclear modification factor decouples from the 
pQCD description. This means that very probably quark coalescence [32-34] dominates 
hadron production in the intermediate momentum region under 6-7 GeV/c, both at RHIC 
and LHC energies, although many questions connected to hadron-hadron correlations 
remained open. 

The new ALICE data at ^/s = 2.76 ATeV [13] induced an extended discussion and 
many papers appeared after the publication of these data (see e.g. [35 40]). Final conclu- 
sion can be drawn after the new pp data collected at ^/s — 2.76 TeV will be analysed and 
direct comparison to data with small statistical and systematical errors can be performed. 



6. Conclusions 

We have investigated jet energy loss at RHIC and LHC energies. Comparing the ex- 
perimental data on nuclear modification factor measured in central Au+Au and Pb+Pb 
collisions, we can claim the presence of a weak nuclear shadowing effect to able to ac- 
commodate a larger jet energy loss at LHC energies with respect to the proper opacity 
values at RHIC energies. In case of weak nuclear modification factor the jet energy loss 
can explain the measured data. These calculations should be repeated when experimental 
data on reference pp-collisions will be available at ^/s = 2.76 TeV and a more precise 
nuclear modification factor will be displayed. On the other hand the decreasing tendency 
of the nuclear modification factor at > 8 — 10 GeV/c clearly displays the presence 
of jet energy loss and these new data at LHC deny the presence of simple participant 
nucleon number scaling, which idea was very popular at RHIC energies to explain the 
very flat nuclear modification factor at high-pr in AuAu collisions . 
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